Thermal background can solve the cosmological moduli problem 
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It is shown that the coherent field oscillation of moduli fields with weak or TeV scale masses 
can dissipate its energy efficiently if they have a derivative coupling to standard bosonic fields in a 
thermal state. This mechanism may provide a new solution to the cosmological moduli problem in 
some special situations [jj. 
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Modern theories of high energy physics contains a number of scalar fields which have a flat potential and interacts 
with ordinary particles only with the gravitational strength In the context of low-energy supersymmetry these 
' moduli fields typically have a flat potential intrinsically and acquire a mass of order of weak or TeV scale when 
supersymmetry is broken 0] , although other mechanisms of moduli stabilization have also been extensively discussed 
f"*^ ' recently Q. 

During inflation in the early universe supersymmetry is spontaneously broken in a different manner than it is 
' today due to the large vacuum energy density. Then moduli fields, which we denote by </>, typically acquire a mass of 
order of the Hubble parameter and they relax to a potential minimum at this stage, which is deviated from today's 
value at the current potential minimum by up to the gravitational scale, A<f> < Mq = 2.4 x 10 18 GeV. After inflation 
their mass is turned off to a much smaller value due to the disappearance of vacuum energy density and they do not 
move until the Hubble parameter decreases to their eventual mass scale which is presumably of order of the weak 
scale or TeV scale as stated above. The scalar fields then start coherent oscillation with the initial amplitude up to 
y—( ' Mg, which will dominate the energy density of the universe eventually. According to the conventional estimates, their 
lifetime is given by 

\o • 

^: r * * 4 « 10 (iWev) sec < w 

o-i. 

that is, they decay after theprimordial nucleosynthesis creating a huge amount of entropy to demolish the successful 
| primordial nucleosynthesis 

In this Letter we present a new class of solution to the above cosmological moduli problem by arguing that the 
previous estimate of the decay rate (P), which has been used in all the other proposed solutions to the problem 
• »"H , does not apply in a finite-temperature and finite-density state of the early universe and that it is much more enhanced 
than in the case of decay in a vacuum. As a result we show that the coherent moduli oscillation may efficiently 
^ i dissipate its energy density well before the big bang nucleosynthesis in some special situations. 

The crucial point is to take moduli decay through derivative coupling correctly into account, such as a coupling with 
a kinetic term of other fields. Indeed we expect moduli field is coupled with gauge fields through -^F^F^ . It may 

also be coupled with scalar fields as \jj^(dx) 2 or A-^x n X where x is a generic scalar field and A is a dimensionless 
coupling constant of order unity Q . Previously, these couplings were expected to give a decay rate no different from 
(fTj) at most, because using the equation of motion, = m x x, it was concluded that the derivative coupling would 
give a decay rate similar to the coupling m^ to xX 2 ' which yields ([TJ for > 2m x . 

Such a naive analysis, however, could only be valid in decay in the vacuum and would not apply in the high- 
temperature environment in the early universe. If x is a standard field in the visible sector, it is strongly coupled with 
other degrees of freedom in the early universe and rapidly reaches thermal equilibrium. Then they acquire a thermal 
mass of order of ~ gT in general, where g is a typical gauge coupling and T is the cosmic temperature, so that the 
right-hand-side of the above equation of motion could be significantly enhanced. Alternatively, one may regard the 
derivative d acting on x as not yielding its rest mass energy m x but energy-momentum arising from finite-temperature 
environment, d ~ T modulo some coupling. Then, the strength of interaction A-p^-(c>x) 2 is estimated as 



(2) 
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The decay rate of <f> through the above interaction should read 

\ 2 g 4 T 4 



where C is a suppression factor due to a large thermal mass of the decay product \ which has been given in 0] for 
a specific model. Here tib(w) = l/(e"/ T — 1) is the thermal number density of a boson and the factor in the bracket 
represents the effect of induced emission 0, 0|. Taking A ~ g ~ 1, ~ 10 2 GeV, T ~ 10 10 GeV, which is a typical 
radiation temperature at the onset of moduli oscillation H ~ m^, we find r ~ 3 x 10 8 C GeV. Thus if C takes an 
appropriate value, <j) can dissipate its energy right after it starts oscillation. 

In order to examine if the above naive estimate is correct, we employ nonequilibrium field theory to calculate the 
dissipation rate of a modulus <f> in the presence of a derivative interaction. For simplicity we consider the following 
model consisting of two scalar fields, <f> and X- 

C = \ (d^) 2 - \ m 2 </> 2 + \ {d, X f - \ m 2 xX 2 + ^{d X f - Ant , (4) 

where Ci nt is interaction of x which thermalizes it. 

We calculate the dissipation rate of (f> under the following setup appropriate to the specific problem we are working 
on. First we neglect cosmic expansion since we are interested only in the case dissipation time is shorter than the 
cosmic expansion time. Second we assume x is in a thermal state with a specific temperature T = f3 owing to the 
rapid thermalizing interaction due to the self coupling. Finally, we consider the situation the parametric resonance 
[lfj | is ineffective which is the case for the modulus mass range of our interest [ll| . 

We calculate an effective action for cf> and derive an equation of motion for its expectation value using the closed 
time-path formalism 12, 13] ■ Although this method has been applied to various cosmological problems by a number of 



authors 1% |15J, to our knowledge, derivative coupling at finite temperature has not been investigated in this context 



yet. The one- loop effective action relevant to dissipation due to the derivative coupling is given by 

d i x(j)A(x)( □ + )(/>c(x) 
d A xd A x'C{x - x')6{t - t')(f> A (x)<j) c (x') + I / d i xd i x'D{x - x'^^x^^x') + ■■■ . (5) 



Here 4> c and </>a are mean and difference of the field variable in the forward time branch (t = —oo to +oo), <f> + , and 
that in the backward time branch (t = +oo to — oo), <f>_, namely, <j) c = (<p + + cf>_)/2 and 0a = <P+ ~ <t>-i respectively. 
<j) + and 0_ should be identified with each other in the end. The kernels in ([5]) are defined by 

C(x - x>) = ^ri^v^Hm [g^Jx - x')G^ 2 (x - x')] , for t - t' > 0, (6) 



D(x-x') = ^ v ^ V ^Re[gF W2 (x-x>)gF V2 (x-x')) 



(7) 



Here G^(x, a;') is a finite-temperature Feynman propagator of field derivatives defined by 

g^(x,x') = (P\Td^ X (x)dMx')\[3)=d x »d x ,„GF(x-x') + i5^5 u0 8 4 (x-x'), (8) 

with G F x {x) = ((3\T X (x)x(x'M. 

The effective action ([5]) is complex-valued as a manifestation of the dissipative nature of the system. We cannot 
obtain any sensible equation of motion by simply differentiating with respect to a field variable because we are dealing 



with a real scalar field and its equation of motion should be real- valued. As shown in 14j , one can obtain a real- valued 
effective action by introducing a random Gaussian variable £(x), which represents fluctuation related to dissipation, 
with the dispersion (^(x)^(x') ) — D{x — x'). As a result the equation of motion for the expectation value </< c (x) is 
given by 



( □ + m 2 ) <j) c (x) + dt' d 3 x'C(x - x')<t> c (x') = £{x) . (9) 
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Hereafter we omit the suffix c. 

As described in 0, Q, this equation can easily be solved using Fourier transform. As a result, we find the dissipation 
rate of zero-mode modulus oscillation is related to the imaginary part of </>'s self energy and given by 



r„i = i 



2mj, 



(10) 



with Co(m,p) being the k = mode of the Fourier transform, C k (uj), of the memory kernel C(x). 

In order to take thermal effects of \ correctly into account, we should use the full dressed propagator, G P u dls ^ [x, x'), 
to calculate C h =a(u) — m ij>) 0- ^ i s obtained by calculating 



l$ Am \x,x') = {(3\Td^)d vX {x')^ v (-i J C mt d A x^ |/3), 



(11) 



using the Matsubara representation [16[ and resummation. As a result, the spatial Fourier mode of the finite- 
temperature dressed propagator is given in terms of the spectral representation as 



G$ dtB) (p,t)=i I ^ ({[l + n B (u)]0(t)+n B (u)0(-t)} 



1 



1 



(w + iT p ) 2 - w 12 (w - iT p ) 2 



(12) 

x + S_r(p) + Tp and T p = — £j(p)/(2cj), where £.r(p) and are real and imaginary 

parts of x's self energy. The real part typically scales as S^(p) ~ g 2 T 2 . 

Inserting (JT2J) into (dHJ), the dissipation rate of the coherent field oscillation is given by 



with rj p = (u),p), oj' 2 = p 2 + in 



X 2 



d 3 p 1 



2 m<j) M 2 J (2tt) 3 4w; 2 



^ \(2p 2 +m 2 T ) 2 [2n B (uj' p ) + l} 



2T r , 



2T r 



(m* - 2uo' p ) 2 + (2T p ) 2 (m + 2u' p ) 2 + (2r p ) 2 



-2(2p 2 + n%fn%{<J p )eHpr p + 4(2p 2 + m 2 T )io' p T p [2n B (Lo' p ) + 1] 



TO0 + 2oj' p 



(m - 2w' ) 2 + (2r p ) 2 (m + 2^i) 2 + (2r p ) 2 



+ 2m A T n B (J p )e p <f3T p 



2m„ 



< + (2Y p ) 2 



(13) 



to the first order in T p . Here = to 2 + is the finite-temperature mass. We find that the first two terms are at 
most of order of T 3 /Mq for T p < T. The last term is also of the same order of magnitude if T p is much larger than 
m<t> 12|; but it may have a different dependence and could be larger for a sufficiently small T p . So we concentrate on 
it hereafter. 

Assuming that there are N degrees of decay modes with the same strength, the last term of (Ti"3"|) reads 



r* = 



iVA 2 



2A/£to0 



dto 
2tt 2 w 



Ju 2 -m 2 minlfujy^P- 
V 1 i 



(2r P ) : 



The last factor is maximal and equal to 1/2 when = 2r p . In this situation (I14p is given by 



NX 2 m 3 T T 



rp4 



M£,m<t, 



(14) 



(15) 



Hence it could be enhanced by a factor of T/m^ compared with other terms. Indeed it can be larger than the Hubble 
parameter, H « T 2 /M G , if 



NX 2 g 



2 n 3 



16tt ~ H ' 



(16) 



where we have replaced mx by gT. Since the right-hand-side is larger than unity in the oscillation regime, from the 
above condition we need N 3> 16tt / (X 2 g 3 ) . That is, assuming the coupling constants are of order unity, N should be 
larger than 50 or so. 

Thus in such a rather special situation, thermal effects may dissipate the moduli oscillation efficiently. For this to 
be the case the imaginary part of the thermal correction to the mass of the decay product T p must be sufficiently 
small to satisfy T p ~ m^. 
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In summary, we have calculated thermal effects on the dissipation of modulus oscillation and have shown that they 
can dissipate its energy efficiently in some special situations. 
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